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We study Rabi oscillations between the bands of an arbitrary biased superlattice in a tight-binding 
model. We reduce the problem to an equation of Whittaker-Hill type and, in absence of any known 
solutions in closed form, discuss different approximations to describe the oscillations between the 
Bloch bands. We identify regimes of weak and strong inter-band coupling and compare predictions 
for these Rabi oscillations to numerical results. 
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I. INTRODUCTION 

The motion of electrons in a periodic potential under 
a constant external force is one of the fundamental phe- 
nomena of solid-state physics and treated in many text 
books Bloch oscillations Q and transitions between 
energy bands via Landau-Zener tunneling [sj are impor- 
tant effects discussed in this context. They have been 
studied in many different contexts, like electrons in semi- 
conductor superlattices [3] , optical waves in photonic lat- 
tices [5|, and Bose-Einstein condensates in optical lat- 
tices [y]. The latter allow a particularly high degree of 
control over various system parameters. 

In the present work, we consider a two-band tight- 
binding model with an additional tilting force and study 
the force-induced coupling between the two Bloch bands. 
This has been an active area of research for quite some 
time, although mainly focused on electronic motion in 
superlattices [7l-[lTj|. For ultra-cold atoms in optical lat- 
tices [H, [l^, nearly closed two-band systems are also 
realisable using different techniques as described, e.g., 
in Besides the possibility of experimental realisa- 

tion, a closed two-band model is also interesting as a 
simple model system. Here, we discuss the two lowest 
energy bands of an optical lattice V{x) = Vb cos(2fcLa;), 
with the wave vector of the optical lattice Ul = 2tt/Xl. 
All energies are measured in units of the recoil energy 
En = fi^fc|/(2?7i) and we set h = 1. This has the ad- 
vantage, that all parameters of the model Hamiltonian 
depend only a single quantity, the depth of the optical 
lattice Vb, and can easily be computed [isl ]. 

The two-band Hamiltonian in the tight-binding ap- 
proximation reads 

+ J2e^\ai){ai\ + e+Mbi\ + V {\ai){bi\ + h.c.) 

(1) 

with = ±A/2 -I- nF. Here, Ta(b) denote the hop- 
ping strength in band a(b), ef^ the on-site energies, and 
V — CqF the local coupling between the energy bands. 
In a realisation with optical lattices, the parameters 
A, Co, Ta, Tf, depend on the depth of the lattice Vb, fig. [l] 
shows this dependence of the Hamiltonian parameters for 
different values of Vb. In the following, we consider the 



external Stark force as a free parameter and our focus 
will lie on the oscillations between the two Bloch bands. 
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FIG. 1: Hamiltonian parameters as a function of the lattice 
depth Vo- Left panel: Average band gap A. Right panel: 
Hopping coefficients Ta (dashed), rj, (dashed-dotted) , and in- 
terband coupling constant Co (solid). 



II. MOMENTUM SPACE ANALYSIS 

We start by changing to Fourier space in two steps: 
first we remove the force term in the time-dependent 
Schrodinger equation by performing a gauge transforma- 
tion and in a second step we will switch to momentum 
space. The applied force breaks the translational invari- 
ance of the system which can be restored by changing 
to the interaction-picture with respect to the external 
force [iBl , such that the Stark term is absent but momen- 
tum p is replaced hy p — Et and the system is explicitly 
time-dependent. We introduce Fourier components 



a{k) = E e'^'a/ & Hk) = E 



(2) 



and obtain from the original non-interacting two-band 
model the following time-periodic two-level Hamilto- 
nian ill 



H{k,t) = 



A 
2 



Ta COs(fc 

V 



Et) V 

I - T6 cos(fc + Et) 



(3) 
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FIG. 2: The two energy bands (solid line) of Eq. in mo- 
mentum space at f = 0. A non-zero force F leads to periodic 
motion within the energy bands. Parameters are for Vb = 4: 
A = 4.39, Ta = 0.062, n = -0.62. 



The spectrum at fixed time is shown in fig. [21 



Transforming a{k, t) 
Ft')dt')] and h{k,t) 



a(A;,t)exp[i(— 
6(fc,t)exp[i(-h^t 



Fa /q cos(fc 
Tfc /g cos(A: 



Ft')dt')], one arrives at a purely ofi^-diagonal Schrodinger 
equation 



idt 



a{k, t) 
b{k,t) 











a{k,t) 
b{k,t) 



(4) 



wheie (f>{k,t) — A-i — Aa;[sin(fc-|-Ft) — sin(fc)] is the phase 
between the two Bloch bands and Aa; = (tq — ti,)/F. 
Where the pure existence of two energy bands allows phe- 
nomena like Rabi oscillations, it is the non-trivial phase 
difi'erence caused by Ax ^ 0, that gives rise to interest- 
ing, non-trivial phenomena for the time-evolution in both 
bands to be discussed below. Please note, that our model 
is slightly different from a tight-binding superlattice, as 
e.g. in [ll|, since we include a local on-site coupling of the 
bands. Let us discuss the special case Tq = t;, first. The 
non-linear time dependence vanishes and we can write 
the system of differential equations as a linear differen- 
tial equation of second order. We differentiate the second 
equation and insert it into the first one, arriving at the 
simple equation a + iAa + V^a = 0. This is easily solved 
and one finds Rabi-likc oscillations 



\a{k,t)\' 



A2 + 4^2 



sin2(A-t/2) 



(5) 



with frequency A — \J A^ -|- AV^ and an amplitude much 
smaller than unity for the parameters chosen here. The 
same procedure can be applied for the general problem, 
but the differential equation to be solved has now time- 
dependent coefficients 



= 5 



A + (tq - Tfc) cos(fc ^-Ft) a H- V'^~a. (6) 



We make an ansatz a = ?/ exp[if2(t)] and choose 
1 



[A + (ja - n) cos(fc + Ft)] , 



(7) 



to remove the second term in eq. ([6]). We finally arrive 
at 



= 2y y 



2V^2 _ _ sin(fc + Ft) 



-3{A+{Ta~n)cos{k + Ft)y 



(8) 



This equation is of the Whittaker-Hill type, described 
for instance in [l7l |. and solutions in closed form are not 
known for the general case. For the rest of this work, we 
will study different approximate solutions of the problem, 
i.e., of eq. ^ and Q. 

A formal solution of eq. Q is given by 



a(fc, t) 
b{k,t) 



— Texp 



-1 / H{k,t')dt' 
'o 



a(fc,0) 
6(fc,0) 



(9) 



where T denotes time-ordering. Following Zhao et al. [8| , 
it can be solved approximately by dropping the time- 
ordering (valid to first order in V). Choosing occupation 
in the lower band only as initial state, i.e. a{k, 0) = 1 
and b{k, 0) — 0, one obtains for the occupation of the 
upper band 



|6(fc,<)p 



V 



(10) 



with 0(fc, i) — A -t — Ax[sin(A; -I- Ft) — sin(fc)] as above. 
Unfortunately, this integral can not be expressed in sim- 
ple terms of elementary functions, but has to be treated 
in different levels of approximation or numerically. To 
obtain an idea of the physics contained in the above 
expression, we stick to some simple approximation and 
will find signs of resonances in this system. Writing 
e"^ = cos <f) + i sin (j) and using the addition theorems, we 
find integrands of the form cos(Ax sin t) and sin(Aa: sin t) 
which are not accessible. We therefore expand the func- 
tions in Fourier series 



cos(A2:sint) = Jo (Ax) -f ^ 2 J2„(Aa::) cos(2nFt) 

n=l 

oo 

sin(A2:sini) = ^ 2J2n+i(Aa;) sin ((2n + l)Ft) 



n=0 



where Ji{x) denotes the ordinary Bessel function of or- 
der I. Taking only the first terms, the integrals can be 
solved analytically and expressions of the following form 
(n = 1,2) appear 



cx sm 



sin"(Aa;/n) 



sin[(A - nF) t] 



AnF 



The higher orders in the expansion lead to similar terms. 
Thus, a truncated Fourier expansion indicates the exis- 
tence of resonances in the occupation of the upper band 
for A « nF, that is for values of the external force where 
the tilt between sites of different bands equals the band 
gap. To test this expectation, we evolved an initial state 
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FIG. 3: Occupation of band 2 as function time in resonance. 
We observe long period and large amplitude Rabi oscillations 
with smaller and much faster oscillations on top. Parameters 
as in fig. [2] and F = 2.2207. 
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FIG. 4: Long time averages of the occupation of band 2 as 
function of the inverse force. There is a clear signature of reso- 
nant behaviour when the band gap equals an integer multiple 
of the force. The dashed line shows the average occupation of 
band 2 as expected from eq. ([5]) and the inset schematically 
shows an n = 2 resonance. Parameters as in fig. (2] 



non-resonant oscillations are well described by eq. (O, 
compare also fig. ID Below we will derive an effective 
model for the system in resonance that allows a quantita- 
tive prediction of the resonance width. Similar resonant 
behaviour has been observed in extensive numerical stud- 
ies of the occupation of minibands for electronic motion 
in superlattices 



III. SYSTEM IN RESONANCE 

Having found resonances in the system, we are now go- 
ing study this resonant behaviour in more detail. Start- 
ing from the original Hamiltonian in coordinate space, 
eq. ([T]), we apply a transformation known [10] to diago- 
nalise the single band problem 



E 



Jl-n{TalF)\ai) 



(11a) 



(lib) 



where we defined Xa(b) '■= Ta(b)/F and Ax — Xa — Xb- 
Using properties of Bessel functions, we arrive at 



n^Y.\ei\ai){ai\+ei\mi\ + 
I 

J,_„(Aa;)(| azXA+nl + h.c.) 



(12) 



We obtain a coupling between any two sites of the two dif- 
ferent bands, weighted by Bessel functions. This means 
the coupling between different and possibly remote sites, 
which is originally mediated by on-site coupling and sub- 
sequent hoppings, can be considered a higher order pro- 
cess in the original basis. But in the transformed basis 
eq. (jl2p . this coupling is now direct with a strength mod- 
ified by the factor Ji_„(Aa;). This is schematically de- 
picted in fig.[S] Interestingly, the same matrix representa- 



with a(fc,0) = 1 and 6(fc, 0) = in time for a resonant 
value and the resulting occupation of the upper band 
|6(i)p is shown in fig. [3l We observe small Rabi oscil- 
lations with a high frequency, which are described by 
eq. ([S]), and on top of these are Rabi oscillations of large 
amplitude and long period T ^ Tg. The latter oscilla- 
tions only take place near resonances A k, nF. 
We did similar time evolutions for various values of the 
external force F and studied the long time average of 
|6(i)p. The result is depicted in fig. ID It clearly demon- 
strates the existence of resonances in our system. The 
occupation of the upper band is usually very small, of 
the order of a few percent, but shows sharp resonances 
for nF « A. A closer numerical study (not shown 
here) shows a Lorentzian shape for the lower order reso- 
nances and the possibility of more complicated, Fano-like 
shapes for the narrower higher order resonances. The 




FIG. 5; Graphical representation of the two-band Hamilto- 
nian: a) in the original basis, eq. ([T]), and b) in the trans- 
formed basis, eq. (|12p . 

tion can be found by application of Floquet theory [l^ to 
the time-periodic Hamiltonian eq. (j?]). One has to solve 
the larger but time-independent eigenvalue problem for 



4 



the Floquet operator H f having matrix elements 



This is based on expanding the Hamihonian in a Fourier 
series with Fourier components . In our case these 
Fourier components are 2x2 matrices and read exphcitly 



H 



2n+l 



-A/2 V 

V A/2 



T/J|2„|(Ax) 



1 

1 



sgn(n)T/J|2„+i|(Ax) 



-1 

1 



(13a) 
(13b) 
, (13c) 



with the signum function sgn(a:). The obtained Flo- 
quet Hamiltonian Hp is the same as obtained earlier in 
eq. (|12p . For the special case Ax — 0, i.e. Ta = t^, the 
model eq. can be solved exactly as before, cf. eq. ([SJ. 
The argument of the Bessel functions is zero for Ax = 
and only on-site coupling of the two bands remains. 

When the system is exactly in resonance, there is a 
degeneracy between two energy levels of the two bands 
and the most relevant coupling in the transformed basis, 
is the direct one between these two levels. Taking only 
this coupling, we arrive at an effective Rabi problem [19| . 
In a resonance of order m, this reduces to a simple 2x2- 
matrix: 



T,(m) _ 



VJm{Ax) e+ 



(14) 



The occupation of the upper band follows a simple Rabi 
formula 



sin^ {n,,,t) , (15) 



(toF- A)2 +4F2J2„ 
with a Rabi frequency for the resonant oscillations given 



by fires = h 



E+)+4F2j2 



|yj„(Ax) 



This frequency is much smaller than the Bloch frequency 
ojb = 27r/rB and the corresponding periods are given by 



F 



2\VJmiAx)\ 



(16) 



Thus, the system in resonance exhibits oscillations be- 
tween the two bands of an amplitude close to unity and 
on a time scale much larger than the direct on-site oscilla- 
tions. In fig. |5]we compare this predictions to numerical 
simulations for different values of the parameters in the 
Hamiltonian, eq. ([Ij, and different orders of resonance. 
We observe very good agreement. Additionally, the inset 
shows the relative error between the predicted and nu- 
merically observed values. It is approximately 5 % for 
the widest resonance n — 1, but smaller than 3 % for 
higher resonances and slightly decreasing with increasing 
potential depth. 
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FIG. 6: Periods of resonant oscillation exactly in resonance 
as predicted by mapping on a two-level Rabi problem (dashed 
lines) and compared to periods from numerical time evolution 
in perturbation theory (data points). The inset shows the 
relative error between the two values. 



We want to compare these results to nearly degenerate 
perturbation theory in the original basis. Here, the res- 
onances of different order have to be treated separately. 
Starting with the first order resonance {A — F), one finds 
in second order degenerate perturbation theory a Rabi 
frequency of 51 = ^\V {ti, — Ta) / A\ reproducing the lead- 
ing order result of ilios"^'' — \V Ji{Ax)\. For the second 
order resonance (A — 2F), one obtains in third order de- 
generate perturbation theory V, = ^\V{t^ + t^)/{8F^)\, 
slightly differing from the leading order behaviour of the 
result in the transformed basis. Higher order resonances 
can be computed in a similar fashion, but the resonant 
basis pT|) is clearly a better starting point for further 
calculations. Indeed, considering the mean occupation of 
the upper band close to a resonance, eq. (jlSp . we find a 
Breit-Wigner shape for the resonances as also seen nu- 
merically 0, namely 



1 



{\bmii/F) = - 



2(m/A-l/F)2+r2„ 



(17) 



where F^ = 2\VJ„riAx)\/{FA) and the width of the res- 
onances is accordingly given by 2r,„. Comparing to nu- 
merical simulations (not shown here) , the effective model 
presented here gives a good estimate for the magnitude of 
the resonance width. In particular, the width decreases 
very quickly with the order of the resonance or the depth 
of the optical lattice, since 



2r 



4|Co| m^' 
A 2™^ 



n 



(18) 



where the factor ni™/2"^ml is of order one (for m — 
1, . . . , 10), but since (r^ — r^) ^ A, the last factor de- 
creases very rapidly. We do not compare these values 
with results from numerical calculations in more detail 
here. Higher accuracy can be obtained by degenerate 
perturbation theory in the resonant basis. 
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IV. SUMMARY 

We studied oscillations in a two-band tight-binding 
model. We could reduce the problem of solving the time- 
dependent Schrodinger equation to a Whittaker-Hill 
equation and concentrated on different physical regimes 
for the inter-band oscillations. We found and discussed 
resonances in the system and derived an effective model 
for the system in resonance allowing a quantitative char- 



acterisation of the resonances themselves and the time- 
scales of inter-band oscillations. 
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